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ABSTRACT 
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O [ In a recent paper (Dai & Lu 1999), we have proposed a simple model in which the 

(N 

o ! steepening in the light curve of the R-band afterglow of the gamma-ray burst (GRB) 

(D ■ 
pL| ' 990123 is caused by the adiabatic shock which has evolved from an ultrarelativistic phase 

P^ ■ to a nonrelativistic phase in a dense medium. We find that such a model is quite consistent 

cn ■ with observations if the medium density is about 3 x 10^ cm~'^. Here we discuss this model 

^ . in more details. In particular, we investigate the effects of synchrotron self absorption 

>Q ! and energy injection. A shock in a dense medium becomes nonrelativistic rapidly after a 

Qv^ I short relativistic phase. The afterglow from the shock at the nonrelativistic stage decays 



o 



more rapidly than at the relativistic stage. Since some models for GRB energy sources 
predict that a strongly magnetic millisecond pulsar may be born during the formation of 
GRB, we discuss the effect of such a pulsar on the evolution of the nonrelativistic shock 
through magnetic dipole radiation. We find that after the energy which the shock obtains 



k> i from the pulsar is much more than the initial energy of the shock, the afterglow decay 

c^ I will flatten significantly. When the pulsar energy input effect disappears, the decay will 

steepen again. These features are in excellent agreement with the afterglows of GRB 

980519, GRB 990510 and GRB 980326. Furthermore, our model fits very well all the 

observational data of GRB 980519 including the last two detections. 

Subject headings: gamma-rays: bursts - stars: pulsars - shock waves 



1. INTRODUCTION 

In the standard afterglow shock model (for a review see Piran 1999), a gamma-ray 
burst (GRB) afterglow is usually believed to be produced by synchrotron radiation or in- 
verse Compton scattering in an ultrarelativistic shock wave expanding in a homogeneous 
medium. As more and more ambient matter is swept up, the shock gradually decelerates 
while the emission from such a shock fades down, dominating at the beginning in X-rays 
and progressively at optical to radio energy band. In this model, there are two limiting 
cases (adiabatic and highly radiative) for hydrodynamical evolution of a relativistic shock. 
These cases have been well studied both analytically (e.g., Meszaros & Rees 1997; Wijers, 
Rees & Meszaros 1997; Waxman 1997a, b; Reichart 1997; Sari 1997; Vietri 1997; Katz & 
Piran 1997; Sari, Piran & Narayan 1998; etc) and numerically (e.g., Panaitescu, Meszaros 
& Rees 1998; Huang et al. 1998; Huang, Dai & Lu 1998). A partially radiative (interme- 
diate) case has been investigated (Chiang & Dermer 1998; Cohen, Piran & Sari 1998; Dai, 
Huang & Lu 1999; Huang, Dai & Lu 1999a). All the studies are based on the following 
basic assumptions: (1) the total energy of the shock is released impulsively before its for- 
mation; (2) the medium swept up by the shock is homogeneous and its density (n) is the 
one of the interstellar medium ~ 1 cm~^; and (3) the electron and magnetic field energy 
fractions of the shocked medium and the index (p) in the accelerated electrons' power-law 
distribution are constant during the whole evolution stage. The standard model is suc- 
cessful at explaining the overall features of late afterglows of some bursts such as GRB 
970508: the light curves behave according to a single unbroken power law with decay 
index of a ~ —1 as long as the observations continued (Zharikov et al. 1998). 

Each of these assumptions has been varied to discuss why some observed afterglows 
deviate from that expected by the standard afterglow model. For example, the R-band 
light curve of GRB 970508 afterglow peaks around two days after the burst, and there is 
a rather rapid rise before the peak which is followed by a long power-law decay. There are 
two models explaining this special feature: (i) Rees and Meszaros (1998) envisioned that 



a postburst fireball may contain shells with a continuous distribution of Lorentz factors. 
As the external forward shock sweeps up ambient matter and decelerates, internal shells 
will catch up with the shock and supply energy into it. A detailed calculation shows that 
this model can explain well this special feature (Panaitescu, Meszaros and Rees 1998). 
(ii) Dai & Lu (1998a) considered continuous energy injection from a strongly magnetized 
millisecond pulsar into the shock through magnetic dipole radiation. This model can 
also account for well the observations. It is very clear that these models don't use basic 
assumption (1). 

There are several models in the literature that discuss the effect of inhomogeneous 
media on afterglows (Dai & Lu 1998b; Meszaros, Rees & Wijers 1998; Vietri 1997), 
dropping the second assumption. Generally, an n oc r~^ {k > 0) medium is expected to 
steepen an afterglow's temporal decay. Recently, Chevalier & Li (1999a, b) found that a 
Wolf-Rayet star wind likely leads to an n oc r~^ medium, and thus if GRB 980519 resulted 
from the explosion of such a massive star, subsequent evolution of a relativistic shock in 
this medium is consistent with the steep decay in the R-band light curve of the afterglow 
from this burst. Another way of dropping the second assumption is that the density of 
an ambient medium is invoked to be as high as n ~ 10^ cm^^. Recent observations show 
that the temporal decay of the R-band afterglow of GRB 990123 steepened about 2.5 
days after this burst (Kulkarni et al. 1999; Castro-Tirado et al. 1999; Fruchter et al. 
1999). Dai & Lu (1999) (hereafter DL99) proposed a plausible model in which a shock 
expanding in a dense medium has evolved from a relativistic phase to a nonrelativistic 
phase. They found that this model fits well the observational data if the medium density 
is about 3x10^ cm~^. They further suggested that such a medium could be a supernova 
or supranova or hypernova ejecta. Of course, the steepening in the light curves of the 
afterglows of these two bursts may be due to lateral spreading of a jet, as analyzed 
by Rhoads (1999) and Sari, Piran & Halpern (1999). However, numerical studies of 
Panaitescu & Meszaros (1999), Moderski, Sikora & Bulik (1999), Huang et al. (1999b) 



and Huang, Dai & Lu (2000) show that the break of the hght curve is weaker and much 
smoother than the one analytically predicted when the light travel effects related to the 
lateral size of the jet and a realistic expression of the lateral expansion speed are taken 
into account. 

In basic assumption (3), the electron and magnetic field energy fractions of the shocked 
medium may not be varied during whole evolution, as argued by Wang, Dai & Lu (1999a), 
who analyzed all the observational data including both the prompt optical flash and 
the afterglow of GRB 990123. However, the assumption that p is constant might be 
inconsistent with the early afterglow from GRB 970508 (Djorgovski et al. 1997). 

In this paper, we discuss the model proposed by DL99 in more details, by taking into 
account both the synchrotron self-absorption effect in the shocked medium and the energy 
injection effect of Dai & Lu (1998a, c). Therefore, our present analysis, in fact, relaxes 
assumptions (1) and (2). So far the bursts whose afterglow decay steepens include GRB 
980519, GRB 980326 and GRB 990510 besides GRB 990123. In particular, for the former 
two of these bursts, optical observations several days later are far above a power law 
decline, implying possible energy injection at such a late stage. In section 2, we analyze 
the spectrum and light curve of radiation from a shock expanding in a dense medium. In 
section 3, we compare our model with observations related to GRB 980519 and infer all 
intrinsic parameters and the redshift of this burst. We discuss properties of GRB 990510 
and GRB 980326 in section 4, and in the final section we give a discussion and conclusion. 

2. SHOCK EVOLUTION 

2.1. Relativistic Phase 

For simplicity, we assume that a relativistic shock expanding in a dense medium is 
adiabatic. The evolution of a partially radiative shock depends on both the efficiency 
with which the shock transfers its bulk kinetic energy to electrons and magnetic fields 
and on the efficiency with which the electrons radiate their energy (Dai, Huang & Lu 



1999). Here we don't consider such a shock. The Blandford-McKee (1976) self similar 
solution gives the Lorentz factor of an adiabatic relativistic shock, 
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where Eq = E^2 x lO^^ergs is the total isotropic energy, n^ = n/10^ cm"^, t^ is the 
observer's time since the gamma-ray trigger in units of 1 day, z is the the redshift of the 
source generating this shock, and nip is the proton mass. We assume 7 = 1 when t© = t&. 
This implies 

n, = E,2t^%l + z)/2f. (2) 

For t® > tb, the shock will be in a nonrelativistic phase. In the following we will see 
different spectra and light curves from this shock before and after the time tb- We first 
analyze the relativistic case. 

As usual, only synchrotron radiation from the shock is considered. To analyze the spec- 
trum and light curve, one needs to know three crucial frequencies: the synchrotron peak 
frequency (z/m), the cooling frequency (uc), and the self-absorption frequency (z/q). We 
assume a power law distribution of the electrons accelerated by the shock: dn'^/d'-fe oc 7^^ 
for 7e > 7em, where 7e is the electron Lorentz factor and '-fem = 610ee7 is the minimum 
Lorentz factor. We further assume that e^ and e^ are the electron and magnetic en- 
ergy density fractions of the shocked medium respectively. The z/^ is the characteristic 
synchrotron frequency of an electron with Lorentz factor of '~fem, while the z/^ is the char- 
acteristic synchrotron frequency of an electron which cools on the dynamical age of the 
shock. According to Sari et al. (1998), therefore, these two frequencies, measured in the 
observer's frame, can be written as 

r^m = yI^^I^ = 7.0 X 10\eJO.ire]i%EU%'^'[{l + z)/2]'/' Hz, (3) 

ISnerUecil + z) ^^ , „i7 -3/2 j^-1/2 -1,-1/2 
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2.2 X lO''eZlE^2^WW^ li^ + ^)/2]-^/' Hz, (4) 



where 6^,-6 = ^s/lO"®, B' = (327reB7^H«^pC^)^^^ is the internal magnetic field strength of 
the shocked medium and a^ is the Thompson scattering cross section. The self-absorption 
frequency has been estimated to be 

z/, = 6.0 X 10'\ee/0.irhf_,Eli'nf[{l + z)/2]-' Hz, (5) 

where p = 2.8 has been used (Granot, Piran & Sari 1998; Wijers & Galama 1998). This 
estimate is valid only for Ua < I'm- Since z/^ decreases as oc t0 and Ua is time invariant, 
we define a time t^ based on Vrn{ia) = ^a'- 

ta = O.OlieJO.l)h]i^_,^Eli%f^'[{l + z)/2] days. (6) 

When Ua > i^m, the self-absorption coefficient a,y oc 77^„^i?'(P+^)/^z/~(P+'')/^ (Rybicki & 
Lightman 1979) and the width of the shock Ar ^ r/'-f oc 7t0, so, based on the optical 
depth r = a^Ar = 0.35, we find that the self-absorption frequency decays as Ua oc 
^^(3p+2)/2(p+4)_ Thus we have 

ua = 6.0 X 10'\ej0.1)-h]i'_Xi'nT[i^ + z)/2]-\tJQ-^^ Hz, (7) 

The observed synchrotron radiation peak flux can be obtained by 

""^ 4^ = 10eB%E52n/ \^—^=——j mJy, (8) 

where N^ is the total number of swept-up electrons and P^^ = meC^(TTB'/{3e) is the 
radiated power per electron per unit frequency in the frame comoving with the shocked 
medium. For a flat universe with Hq = 65kms~^ Mpc~^, the distance to the source 
Dl = 2c/Ho{1 + z-^/TT^). 

After having the peak flux and three break frequencies, we can write the spectrum and 
light curve of synchrotron radiation. For high frequency u > Vam = niax(z/(j, z/^), we flnd 

p ^ / (^/^™)-(^-^)/^F.. oc z/-(P-^)/2tf -^)/^ if z/_ < z/ < z/, 

"" I {uJUm)-^^-^y\u/u:)-^I^F^^ oc U-^lHf^" if Z. > Z/e. ^^ 



If p ^ 2.8, then the temporal decay index a = 3(1 — p)/4 ?a —1.35 for emission from 
slow-coohng electrons or a = (2 — 3p)/4 ^ —1.6 for emission from fast-cooling electrons. 
In addition, the low-frequency (z/ < Uam) radiation should be discussed in two cases: (i) 
for Ua < I'm, the spectrum and light curve can be written 

I {vjVmYl'^F^^ OC z/^/^tJ/^ \iva<v < z/„; 

(ii) for Va > Vmi we can obtain the spectrum and light curve, 

F. = {vJv^„:)'^^-^^l\vlv:fl^F^^ OC v^lH't. (11) 

These equations show that the flux of the low-frequency radiation increases with time. 

2.2. Nonrelativistic Phase 

After sweeping up sufficient ambient matter, the shock will eventually go into a non- 
relativistic phase, viz., t® > 1^. In the following we analyze the spectrum and light curve 
of the synchrotron radiation from such a shock, by assuming z/^ > Um- 

2.2.1. Without Any Energy Injection 

We ffist consider the widely-studied case without any energy to be input into the shock 
after the GRB. In this case, the shock's velocity f oc t^ and its radius r ^ ft^ oc t^ . 
According to DL99, thus, the synchrotron peak frequency, the cooling frequency, the 
self-absorption frequency, and the peak flux are derived as 

Urn = 7.0 X 10\ee/0.lfe]i'_,EU%'/\l + z)/2f'{tjt,r' Hz, (12) 
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2.2 X lO''e-'LlE^,'^\^%'^'[il + z)/2r'/'it^/t,)-'/' Hz, (13) 



Ua = 6.0xlO'\ej0.1)-hf_^EU''nl/''[{l + z)/2]-\tb/tar^^{t^/tbr^ Hz,(14) 



and 

F.^ = 104/%E5.nf (I)''' [ ^'_\' ) " -Jy- (15) 

Based on these equations, we further derive the spectrum and hght curve, 

{u/u^)'(^'^y^F,^ oc ^-(P-1)/2421-I5p)/10 if ^^ < ^ < ^^ (16) 

We easily see that for high-frequency radiation the temporal decay index a = (21 — 15p)/10 
for emission from slow-cooling electrons or a = (4 — 3p)/2 for emission from fast-cooling 
electrons, li p ^ 2.8, then a ~ —2.1 or —2.2. Comparing this with the relativistic result, 
we conclude that the afterglow decay steepens at the nonrelativistic stage. 

2.2.2. With Energy Injection from Pulsars 

Some models for GRB energy sources (for a brief review see Dai & Lu 1998c) predict 
that during the formation of an ultrarelativistic fireball required by GRB, a strongly 
magnetized millisecond pulsar will be born. If so, the pulsar will continuously input its 
rotational energy into the forward shock of the postburst fireball through magnetic dipole 
radiation because electromagnetic waves radiated by the pulsar will be absorbed in the 
shocked medium (Dai & Lu 1998a, c). Since an initially ultrarelativistic shock discussed 
in this paper rapidly becomes nonrelativistic in a dense medium, we next investigate the 
evolution of a nonrelativistic adiabatic shock with energy injection from a pulsar. The 
total energy of the shock is the sum of the initial energy and the energy which the shock 
has obtained from the pulsar: 

Eq+ r Ldt^ = Etot oc vh-^, (17) 

Jo 

where L is the stellar spindown power oc (1 + te/T)"^ (T is the initial spindown time 
scale). The term on the right-hand side is consistent with the Sedov solution. Please 



note that L can be thought of as a constant for t^ < T, while L decays as oc tj^^ for 
t® ^ T. Because of this feature, we easily integrate the second term on the left-hand side 
of equation (17). We now define a time at which the shock has obtained energy ~ Eq 
from the pulsar, tc = Eq/L, and assume tc <^ T. We next analyze the evolution of the 
afterglow from such a shock at three stages. 

First, at the initial stage, t^ <^ t^, viz., the second term on the left-hand side of 
equation (17) can be neglected. The evolution of the afterglow is the same as in the 
above case without any energy injection. 

Second, for T > t^ ^ tc, the term Eq in equation (17) can be neglected. At this stage, 
the shock's velocity f oc t^ , its radius r oc t^ , the internal field strength B' cc t^ 
and the typical electron Lorentz factor 7em oc t^ . Thus, we obtain the synchrotron 
peak frequency Vm oc 7em-B' oc t^, the cooling frequency Vc oc B'^^t^ oc t0 , the self- 
absorption frequency Va oc t® , and the peak flux Fy^ oc N^Pl^ oc r^B' oc t^' . 
According to these scaling laws, we derive the spectrum and light curve of the afterglow 

(z//z/„)-(P-l)/2F.^ OC ^-(P-l)/2t£2-5p)/5 if ^^ < ^ < ^^ (18) 

(^c/^.0-^^-'^/'(^/^c)-^/'F.™ OC V-^/Hl^ if Z/ > Vc 

It can be seen that for high-frequency radiation the temporal decay index a = (12 — 
5p)/5 ~ —0.4 for emission from slow-cooling electrons or a = 2 — p ^ —0.8 for emis- 
sion from fast-cooling electrons if p ^ 2.8. This shows that the afterglow decay may 
significantly flatten due to the effect of the pulsar. 

Third, for t^ ^ T, the power of the pulsar due to magnetic dipole radiation rapidly 
decreases as L oc t^"^, and the evolution of the shock is hardly affected by the stellar 
radiation. Thus, the evolution of the afterglow at this stage will be the same as in the 
above case without any energy injection. 

3. OBSERVED AND INFERRED PARAMETERS OF GRB 980519 



We have shown that as an adiabatic shock expands in a dense medium from an ultra- 
relativistic phase to a nonrelativistic phase, the decay of radiation from such a shock will 
steepen, subsequently may flatten if a strongly magnetic millisecond pulsar continuously 
inputs its rotational energy into the shock through magnetic dipole radiation, and finally 
the decay will steepen again due to disappearance of the stellar effect. We next show that 
these effects can fit very well the observed afterglow of GRB 980519. 

GRB 980519 was one of the brightest of the bursts detected by the BeppoSAX satellite 
(Muller et al. 1998; in 't Zand et al. 1999). The BATSE measured fluence above 25 
keV was (2.54 ± 0.41) x 10~^ergscm"^, which places it among the top 12% of BATSE 
bursts (Connaughton 1998). An X-ray afterglow was detected by the BeppoSAX NFI 
(Nicastro et al. 1999). The optical afterglow ~ 8.5 hours after the burst presented the 
most rapid fading of the well-detected GRB afterglows except for GRB 990510, consistent 
with t^^-^^^^-^"^ in BVRI (Halpern et al. 1999), while the power-law decay index of the 
X-ray afterglow, ax = 2.07 ±0.11 (Owens et al. 1998), in agreement with the optical. 
The spectrum in optical band alone is well fitted by a power low //-i-20±o.25^ while the 
optical and X-ray spectra together can also be fitted by a single power law of the form 
^-i.o5±o.io^ In addition, the radio afterglow of this burst was observed by the VLA at 8.3 
GHz, and its temporal evolution oc t%^^^-^ between 1998 May 19.8UT and 22.3UT (Frail, 
Taylor & Kulkarni 1998). 

We now analyze the observed afterglow data of GRB 980519 based on our model. We 
assume that for this burst, the forward shock evolved from an ultrarelativistic phase to 
a nonrelativistic phase in a dense medium at ~ 8 hr after the burst. So, the detected 
afterglow, in fact, was the radiation from a nonrelativistic shock. This implies 7 ~ 1 at 
tft ~ 1/3 days. From equation (2), therefore, we find 

n5~27^52[(l + ^)/2]^ (19) 

If p ^ 2.8, and if the observed optical afterglow was emitted by slow-cooling electrons and 
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the X-ray afterglow from fast-cooling electrons, then according to equation (16), the decay 
index a^ = (21 — 15p)/10 ~ —2.1 and ax = (4 — 3p)/2 ^ —2.2, in excellent agreement 
with observations. Furthermore, the model spectral index at the optical to X-ray band 
and the decay index at the radio band, (3 = —{p — l)/2 ^ —0.9 and a = 1.1, are quite 
consistent with the observed ones, —1.05 ± 0.10 and 0.9 ± 0.3, respectively. 

We next continue to take into account three observational results. First, on May 
21.6UT, the Keck II 10m telescope detected the R-band magnitude R = 23.03 ± 0.13, 
corresponding to the flux Fr ~ 3.5;uJy at t® ^ 2 days (Gal et al. 1998). Considering this 
result in the second sub-equation of (16) together with equations (12), (15) and (19), we 
can derive 

1.8 n.95 pl.95 fl+Z\'-'' ( VTT~Z-1 \-' 

where p ~ 2.8 has been assumed. 

Second, the BeppoSAX observed the X-ray (2-10 keV) flux Fx ~ 1.3 x lO^^yuJy at 
t® ~ 0.65 days (Nicastro et al. 1999). Inserting this result into the third sub-equation of 
(16) together with equations (12), (13), (15) and (19), we can also derive 



-2 



Third, the VLA detected the radio flux Fs.sghz ~ 102±19/iJy on May 22.3UT in 1998 
(Frail et al. 1998). This may result from the self-absorption effect in the shocked medium. 
Thus, combining it with the flrst sub-equation of (16) together with equations (12), (14), 
(15) and (19) leads to 
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ere-s^'E^i-^' C-^y''^' f ^ . M -3.2. (22) 



V2-1 
In addition, the total energy of the adiabatic shock, Eq, is approximately equal to the 
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one released initially in gamma-rays (Piran 1999). This implies 

From equations (19)-(23), we infer intrinsic parameters of the shock and the redshift 

of the burst: 

ee~0.16, eB,-6~280, 

^52-0.27, n5~3.4, (24) 

z ~ 0.55. 

Our inferred value of Cg is near the equipartition value, in agreement with the result of 
Wijers & Galama (1998) and Granot, Piran & Sari (1998), while our e^ is also close to 
the value inferred from GRB 971214 and GRB 990123 (Wijers & Galama 1998; Galama 
et al. 1999; Wang et al. 1999a). After considering these reasonable parameters, Wang, 
Dai & Lu (1999b) numerically studied the trans-relativistic evolution of the shock and 
found that our dense medium model can provide an excellent fit to all the observational 
data of the radio afterglow from GRB 980519 shown in Frail et al. (1999). 

If the late afterglow of GRB 980519 had still decayed according to equation (16), the 
inferred R-band fiuxes on the 60th day and 66th day would have been nearly two orders 
of magnitude smaller than the observed values. This would lead to the argument that 
the emission on these two days came from the host galaxy of the burst (Sokolov et al. 
1998; Bloom et al. 1998). We note that, despite of excellent seeing conditions on the 
Keck II telescope. Bloom et al. (1998) found little evidence for extension expected of a 
host galaxy. This implies that there may exist some mechanism by which the shock at the 
late stage had been renewed. As suggested in the above section, this mechanism is that 
a strongly magnetized millisecond pulsar had supplied its rotational energy to the shock 
through magnetic dipole radiation. We can see from the second sub-equation of (18) that 
when tc <^ t® < T, the R-band afterglow decay index a^ = (12 — 5p)/5 ~ —0.4, where 
p ^ 2.8 has been assumed. Combining this result with the observed fiux F^ ~ 0.2 /iJy 
on the 60th day and the decay power law in several days after the burst, we infer tc ~ 4 
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days. According to the definition of tc = 29-E52-B7i3-Pms days where -Bs,i3 is the surface 
magnetic field strength of the pulsar in units of 10^^ G and P^s is its initial period in 
units of 1 ms (Dai & Lu 1998a, c), we can obtain a constraint on the stellar parameters: 
Bs,n ~ 2.7-^52 -Pjns ~ ^-''^-^ms- Morcovcr, our model requires T > 66 days, which leads 
to Pjns < 0.8, where we have used the definition of the stellar spindown timescale, T = 
120B~f^P^^ days. Therefore, if GRB 980519 resulted from a pulsar, and if the property of 
the late afterglow was caused by the effect of the stellar magnetic dipole radiation, then 
this pulsar may be a strongly magnetized millisecond or even submillisecond one. 

4. PROPERTIES OF OTHER BURSTS 

4.1. GRB 990510 

GRB 990510 was detected by the BeppoSAX Gamma-Ray Burst Monitor (Piro et al. 
1999) as a bright and complex GRB composed by two well seperated and multi-peaked 
pulses with a total duration of about 75 s (Amati et al. 1999). Its fluence was among 
the highest of the BeppoSAX localized events, after GRB 990123, GRB 980329 and GRB 
970111. It was also detected by BATSE (Kippen et al. 1999) and its fiuence (> 20 keV) 
was (2.56 ±0.09) x lO^^ergcm"^, in the top 9% of the BATSE burst fiuence distribution. 
The burst appears aX z > 1.62 (Vreeswijk et al. 1999), which leads to an isotropic energy 
of > 1.4 X 10^^ ergs. The burst's afterglow was detected and monitored at X-ray and 
optical bands. Even though the X-ray afterglow decay light curve is not unlike that seen 
previous X-ray afterglow decays (Kuulkers et al. 1999), the optical afterglow displays its 
special feature: all the temporal decays at V7?/ bands steepened about 1.2 days after the 
burst (Harrison et al. 1999; Stanek et al. 1999; Bloom et al. 1999a; Marconi et al. 1999). 
Initially the optical decay index ai = —0.82 ± 0.02, but about 1.2 days later the index 
became 0:2 = —2.18 ± 0.05. The consistency of ai, 0:2 and the breaking time means that 
the breaking is wide band. This is the first clear observation of a wide band break (Bloom 
et al. 1999a; Harrison et al. 1999). 
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One simple interpretation for this steepening seems that we have been seeing evidence 
for a spreading jet (Rhoads 1999; Sari et al. 1999; Bloom et al. 1999a). As shown 
numerically in Panaitescu & Meszaros (1999), Moderski et al. (1999) and Huang et 
al. (1999b, 2000), however, the evolution of a spreading jet may not lead to a marked 
steepening. Another possible interpretation is that the effect of a strongly magetized 
millisecond pulsar on the evolution of a nonrelativistic adiabatic shock in a dense medium 
has been becoming unimportant. If initially the pulsar was able to change the evolution of 
the shock for GRB 990510 through magnetic dipole radiation, and if the optical afterglow 
came from fast-cooling electrons in the shocked medium, then according to the second 
sub-equation of (18), the temporal decay index a = 2 — p ^ —0.82 ± 0.02. This requires 
p ^ 2.8, which is quite consistent with the value inferred from GRB 980519. When 
the effect of the pulsar on the shock disappeared, the optical afterglow decayed based 
on the third sub-equation of (16), viz., a = (4 — 3p)/2 ^ —2.2, in excellent agreement 
with the observations. Furthermore, the observed breaking time should be equal to the 
stellar spindown timescale in our model, which constrains the pulsar's field strength: 
Bs,i3 ~ lOPms- Thus, the central engine of this burst could be a millisecond magnetar. 

4.2. GRB 980326 

The afterglow of GRB 980326 also had a rapid decline. Groot et al. (1998) derived a 
temporal decay index of a = —2.1 ± 0.13 and a spectral index of /5 = —0.66 ± 0.7 in the 
optical band. This initial decay index, which is similar to that of GRB 980519, suggests 
the evolution of a nonrelativistic adiabatic shock in a dense medium. There is another 
observational result similar to the case of GRB 980519: the decay of the observed optical 
afterglow began to fiatten about 5 days after the burst; this is not the contribution 
of the host galaxy because it is not present at a later time (Bloom & Kulkarni 1998; 
Bloom et al. 1999b). Consequently, the late afterglow might be interpreted as a different 
phenomenon. Bloom et al. (1999b) suggested that this late afterglow could result from 
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a supernova associating with GRB 980326. In our model, this is understood to be the 
emission from the nonrelativistic shock to which has been input energy by a strongly 
magnetized millisecond pulsar. According to the light curve shown in Figure 2 of Bloom 
et al. (1999b), we infer tc ~ 5 days. Thus, the surface field strength of the pulsar could 
be Bs,i3 ~ 2AEli^Pl^. 

5. DISCUSSIONS AND CONCLUSIONS 

There have been two kinds of plausible models for GRB energy sources in the litera- 
ture: one relating to pulsars and another to stellar- mass black holes. Dai & Lu (1998c) 
have summarized possible progenitors involving strongly magnetic millisecond pulsars: 
accretion-induced collapses of magnetized white dwarfs, mergers of two neutron stars if 
the equation of state for neutron matter is moderately stiff to stiff, and phase transitions 
of neutron stars to strange stars. The pulsar progenitor models also include the birth of 
magnetars in supernova explosions suggested by Wheeler et al. (1999). The rotational 
energy of such pulsars at the moment of their formation is as high as a few x 10^^ ergs. 
The efficiency of transformation of the rotational energy to the energy of a relativistic 
outflow and then to the energy of high-frequency radiation may be as high as almost 100% 
(Usov 1994; Blackman, Yi & Field 1996). Such pulsars have been suggested to generate 
possibly anisotropic outflows (Dai & Lu 1998a; Smolsky & Usov 1999), and thus may 
explain the energetics of GRBs, including GRB 990123 as an extreme event if the energy 
flux from the source at the line of sight is only about ten times more than the energy 
flux averaged over all directions. Two by-products in this kind of models are millisecond 
pulsars and relativistic forward shocks forming during the collision of the outflows with 
ambient media. It is natural to expect that the central pulsars affect the evolution of 
postburst shocks and in turn the afterglows from these shocks through magnetic dipole 
radiation. 

In the second kind of GRB source models. Fryer, Woosley & Hartmann (1999) have 
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summarized possible progenitors involving black hole accretion disks: neutron star-neutron 
star binary mergers, black hole-neutron star mergers, black hole-white dwarf mergers, mas- 
sive star core collapses, and black hole-helium star mergers. This kind of models should 
include the supranovae proposed by Vietri & Stella (1998). 

Some of the source models mentioned above, e.g., failed supernovae (Woosley 1993), 
hypernovae (Paczyhski 1998), supranovae (Vietri & Stella 1998), phase transitions of 
neutron stars to strange stars (Dai & Lu 1998a; DL99) and birth of magnetars (Wheeler 
et al. 1999), may lead to dense media prior to the occurrence of GRBs. In addition, 
dense media have also been discussed in the context of GRBs. For example, Katz (1994) 
suggested collisions of relativistic nucleons with a dense cloud as an explanation of the 
delayed hard photons from GRB 940217. DL99 discussed the evolution of an adiabatic 
shock in a dense medium to explain the steepening feature in the light curve of the R-band 
afterglow from GRB 990123. 

Based on these arguments, following DL99, we discuss the evolution of an adiabatic 
shock expanding in a dense medium from an ultrarelativistic phase to a nonrelativistic 
phase in more details in this paper. In particular, we discuss the effects of synchrotron 
self absorption and energy injection on the afterglow from this shock. In a dense medium, 
the shock becomes nonrelativistic rapidly after a short relativistic phase. This transition 
time varies from several hours to a few days when the medium density is from 10^ to a few 
xlO^cm"^, and the shock energy from 10^^ to 10^^ ergs. The afterglow from the shock 
at the nonrelativistic stage decays more rapidly than at the relativistic stage, while the 
decay index varies from —1.35 to —2.1 if the spectral index of the accelerated electron 
distribution, p = 2.8, and the radiation comes from those slow-cooling electrons. Since 
some models mentioned above predict that a strongly magnetic millisecond pulsar may 
be born during the formation of GRB, we also discuss the effect of such a pulsar on the 
evolution of a nonrelativistic shock through magnetic dipole radiation, in contrast to the 
case discussed in Dai & Lu (1998a, c). We find that after the energy which the shock 
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obtains from the pulsar is much more than the initial energy of the shock, the afterglow 
decay will flatten significantly and the decay index will become —0.4. When the pulsar 
energy input effect disappears, the index will still be —2.1. These features are in excellent 
agreement with the afterglows of GRB 980519 and GRB 980326. Furthermore, our model 
fits very well all the observational data of GRB 980519 including the last two detections. 
Of course, if an afterglow of our interest comes from fast-cooling electrons in the shocked 
medium, the decay index of this afterglow will first be —1.6 during the relativistic phase 
and subsequently —2.2 at the nonrelativistic stage in the case oi p = 2.8. If the energy 
input effect of the pulsar during this stage becomes very important, the decay index will 
be —0.8. When this effect disappears, the index will become —2.2 again. The latter values 
of the index are quite consistent with the observations of GRB 990510. This requires that 
the time before the pulsar is able to affect the evolution of the shock is only as short as a 
few hours. It should be pointed out that whether and when a pulsar significantly affects 
the evolution of a shock largely depends upon the following three parameters: the shock's 
initial energy, the dipole magnetic field and period of the pulsar. 

The fiattening of the late optical afterglow light curves of GRB 980519 and GRB 
980326 has also been interpreted as being due to Ib/c supernovae, which is based on the 
obvious reddening of the observed afterglow spectrum. In our model, this fiattening is 
due to the energy input effect of millisecond pulsars through magnetic dipole radiation, 
and the expected afterglow spectrum would be the typical synchrotron spectrum without 
any dust effect. If dust exists in the vicinity of the pulsars, however, dust sublimation 
may lead to reddening of the afterglow spectrum as suggested by Waxman & Draine 
(1999). In explaining the steepening of the afterglow light curve of GRB 990510, we 
have envisioned the disappearance of the pulsar energy input effect. In our model, radio 
afterglows first brighten as oc t^ with the pulsar energy input effect, and subsequently 
fade down as oc t^ ^ when this effect becomes unimportant and the observed frequency 
is smaller than the synchrotron self-absorption frequency. Such a steepening has also been 
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understood to be due to the lateral expansion of relativistic jets (Rhoads 1999; Sari et al. 

-I jn 

1999). Thus, radio afterglows from jets first decay as oc t0 and then as oc f^ (Sari et 
al. 1999). Therefore, there is an obvious difference in radio afterglow light curves between 
the pulsar energy input model and the jet model. We expect that future observations will 
provide evidence for or against the pulsar energy input model. 

In summary, following DL99, we propose a model for several afterglows, in which a 
shock expanding in a dense medium evolves if its central engine is a strongly magnetized 
millisecond pulsar. We show that this model explains well the features of the afterglows 
from GRBs 980519, 980326 and 990510. 
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